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Abstract 

We report on measurements of the total cross section for e + e~ — > 77 for center-of- 
mass energies between 57.4 and 59.5 GeV, using the AMY detector at the TRISTAN 
collider. We set new limits on the production of a possible new s-channel resonance 
decaying into photon pairs. 



1. Introduction 

Recently, the L3 Collaboration at LEP reported on an apparent excess of e + e~ — > 
£ + £~77 events (£ = ji or e) in which the invariant mass of the photon pair was 
clustered around 59 GeV. 111 The L3 Collaboration did not speculate on the origin of 
their events beyond noting that the probability of observing such clustering in any 
5 GeV wide mass bin above 40 GeV was (9(10~ 3 ), if the events were due purely to 
QED. The three other experiments at LEP have also reported on such events. Out 
of fifteen e + e~ — > £+£^77 events found at LEP with M 77 > 40 GeV, five events were 
found with a 77 invariant mass in a 1 GeV range near 59 GeV. 

A few speculative models have been advanced 13 ' 41 on the origin of the clustered 
events of L3 assuming that they are not due to QED. The models are based on the 
decay of the Z° into a massive object X and another boson, where X may be a scalar- 
pseudoscalar mixture or a spin-2 particle with a mass of approximately 60 GeV that 
decays predominantly into two photons. Regardless of the details of such models, if 
the particle X can couple with electrons, then the direct s-channel production may 
be observable in the TRISTAN energy range. If there is no coupling of X to fermions, 
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it can still be produced via an s-channel photon that results in X — > 77 plus a soft 
monochromatic photon. Previous studies of e + e~ — > 77,777 at TRISTAN lJ ' 61 were 
consistent with the QED prediction, however, no direct search for a narrow resonance 
was conducted. 

We have searched for the direct production of a new state X via the reaction 
e+e" -> X -> 77, using the AMY detector at the TRISTAN e+e" collider, for 
center-of-mass energies yfs between 57.4 and 59.5 GeV in 250 MeV steps, with 1 to 
2 pb _1 of integrated luminosity per point (19 pb- 1 at 57.8 GeV). 

2. The AMY Detector 

The AMY detector is a general purpose solenoidal-type instrument employing 
two inner tracking chambers (VTX and ITC), a central drift chamber (CDC), and 
an electromagnetic calorimeter (SHC), all contained within a 3 Tesla magnetic field, 
and a barrel muon detector (MUO) outside the magnet return yoke. The end cap re- 
gions are instrumented with calorimeters (ESC) and small-angle luminosity monitors 
(SAC). The AMY detector has been described in detail elsewhere. 

Final states with two or more photons are detected primarily by the SHC. This 
is a 14.4 radiation length gas-sampling calorimeter, consisting of twenty alternating 
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layers of lead and resistive plastic proportional tubes. Each tube layer is sandwiched 
between cathode strips in the and 9 directions. The SHC is constructed in six 
sextants, each covering an azimuthal range of 60° and a polar angle range of | cos#| < 
0.73. 

Measurement of electromagnetic shower locations and energies in the SHC is de- 
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termined by the signals from the cathode strips. The strips are arranged in projective 
towers that subtend an angle of 14.2 mr from the interaction point, and are ganged 
internally into five longitudinal divisions. The empirical energy resolution is found 
to be AE/E = 0.23/ ' \[E + 0.06, with E in GeV. The large constant term is thought 
to be due to the effects of operating in the 3 T magnetic field. Shower centroids are 
typically located in and 9 to less than one strip width. 

Triggering 

Within a given layer of the SHC, anode signals are ganged into 48 azimuthal 
sections; the ganged signals are summed externally over multiple layers to form trigger 
signals for 48 towers with either no segmentation or four longitudinal divisions. Three 
redundant combinations of these signals provide triggering for e + e~ — > 77 events: (1) 
total energy trigger — the analog sum of all 48 trigger signals exceeding 8 GeV, (2) 
two-cluster trigger — two separate towers in coincidence with energies exceeding 5 GeV 
for one and 3 GeV for the other, and (3) shower shape trigger — coincidence between 
any two successive longitudinal divisions in a given tower, or in neighboring towers, 
with total energy exceeding 3 GeV. 

Bhabha electrons are used to check the triggering efficiency independently, since 
they behave similarly to 77 events in the SHC. Bhabha events are triggered inde- 
pendently by the SHC triggers and by charged-track triggers derived from signals 
in the CDC and ITC. By examining the frequency of Bhabha events with a track 
trigger but no SHC trigger, we found the SHC trigger efficiency for Bhabha events to 
be essentially 100%. We assume, therefore, that the SHC triggering efficiency for 77 
events exceeded 99% for all of the data of this study. 

4 



Luminosity Determination 

To search for an anomaly in the production of 77 events, the most critical pa- 
rameter aside from event statistics is the reliable determination of the integrated 
luminosity. This is accomplished in the AMY experiment by recording small-angle 
Bhabha events in the ESC, a sampling calorimeter constructed of alternating lay- 
ers of lead and scintillator, with proportional tubes at the depth of electromagnetic 
shower maximum for position determination. The ESC covers the polar angle range 
0.799 < I cos#| < 0.982. The systematic error in ESC luminosity is estimated to be 
~ 2%, and is dominated by the precision on alignment and fiducial definition. 

3. Event Selection and Monte Carlo 

The selection of e + e~ — > wy events (n > 2) was performed according to the 
following criteria: (1) at least two SHC clusters with energies greater than E(, earn /3 
each, (2) the polar angle for each such cluster in the range 45° < 9 < 135°, (3) the 
acollinearity angle between the two most energetic clusters of at most 10°, and (4) no 
charged tracks in the event. The resulting sample of 1054 events was visually scanned, 
resulting in the rejection of 19 background events from cosmic rays, misidentified 
Bhabhas, and SHC noise. 

The cross section for e + e~ — > 77, including radiative corrections to 0(a 3 ), was 
calculated by Monte Carlo integration using a program by Fujimoto, Igarashi, and 
Shimizu! 9 ' 10 ' This program is an implementation of the BASES/SPRING Monte Carlo 
package by Kawabata!" 1 The cross section of (34.81 ± 0.18) pb is based on a fully 
simulated sample of 1.06 fb -1 at 58 GeV that was subjected to the same selection 
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criteria as the experimental data. This value includes an event selection efficiency of 
5 77 = 0.933 ±0.011. 

We estimate the error in normalization for this analysis to be 2.3%, taking into 
account the uncertainty in the luminosity measurement and including the error in the 
77 event selection efficiency. 

4. Results 

Table 1 summarizes the results of the energy scan. The center-of-mass energies 
represent their actual values at the AMY interaction point rather than the nominal 
accelerator energies. The quantity a E represents the rms width of the center-of- 
mass energy distribution. The luminosities L are those obtained with the ESC. The 
quantity iV 77 [Nqed] represents the number of observed [simulated] ny (n > 2) 
events that pass the event selection criteria; Nqed is derived from the Monte Carlo 
simulation at 58 GeV discussed earlier, scaled by s and adjusted to the tabulated 
luminosity. The ratio of the observed 77 events to that expected from the Standard 
Model is also shown for each center-of-mass energy. 

Barrel Bhabha events obtained concurrently with 77 events are used to test the 
internal consistency of the runs and to provide an additional check for systematic 
errors. For each center-of-mass energy, the ratio of observed barrel Bhabha events 
to ESC luminosity is scaled by s. These quantities are given in the last column of 
Table 1, and indicate the constancy in the operation of the SHC during the scan runs. 
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5. Analysis 



The ratio of observed wy events (n > 2) to the QED expectation is plotted 
as a function of center-of-mass energy in Figure 1. The vertical error bars on the 
ratios are statistical only; the horizontal error bars indicate the rms spread in ^/s 
rather than the error on the central values. A comparison of the data to QED alone, 
allowing the normalization to float without constraint, gives a fully consistent fit with 
X 2 /dof = 0.68 for 8 degrees of freedom and a normalization value of 0.961 ± 0.029. 
This fit is shown as the solid curve in Figure 1. The deviation of the normalization 
from unity is consistent, given the 2.9% statistical uncertainty from the fit on this 
parameter as well as the 2.3% normalization uncertainty mentioned above. 

To examine the effect of a new state X, we again allow the normalization of the 
77 data to float without constraint, and then compare these data to a model in which 
the 77 events are produced by an s-channel resonance X of mass M and total width 
T that sums incoherently with the conventional QED processes. 

The effective cross section for production of X at center-of-mass energy y/s, in- 
tegrated over the solid angle of the SHC, is parameterized by a Breit-Wigner form: 

a x (s) = (2 J + l)r ee BR 77 /(s, M, T) , (1) 

where 

/(s,M,r) = (tte ). ff ; (M _^ 2 + r2/4 . (2) 

r ee is the partial width for X — > e + e~, BR 77 = r 77 /T is the branching ratio for 
X — > 77, and (ffc) e ff is the effective acceptance for detecting e + e~ — > 77. We use 
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(fie)eff = ^shc " £77 — 0.707 • 0.933 = 0.659. (This must be modified somewhat if 
the differential cross section is not isotropic.) The cross section a x is then convolved 
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with the radiatively corrected beam resolution function 



t 



G R (s, E) = (^j -jL- j x'- 1 exp [-(z - xf/2] dx , (3) 

E 

where z = (y/s — E)/a E , t = 2(ai/7r)[m(s/m e ! ) — 1], and o E is the width of the 
center-of-mass energy distribution from Table 1, to give the observed effective cross 
section 

oo 

a x (s) = (2 J + l)r ee BR 77 J f(E 2 , M, V) G R (s, E) dE ■ (4) 

o 

The number of events expected from this process at center-of-mass energy for 
an integrated luminosity L is given by Nx(s) = & x (s) ■ L, and adds to the number 
expected from QED alone, Nqed{s), listed in Table 1. 

For a normalization A , the likelihood of observing Nx(si) events from s-channel 
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production of X at center-of-mass energy JJl is a rescaled Poisson distribution: 



C(si) = A exp {Nqed{si) ~ A [N x {si) + N QE D{si)}} 

A [Nx( Sl ) + N QED ( Si )]\ N rr( s i) (5) 



Nqed(si) 

and the overall likelihood is C — Yii £(si) ■ 

This likelihood function peaks at M = (58.35±g;g|) GeV, T = (16 ± 13) MeV and 
(2 J + l)r ce BR 77 = (4.2±|;i) keV, with a normalization of A = 0.955 ± 0.040. The 
best-fit curve of A (Nx + Nqed)/Nqed is shown as the dashed curve in Figure 1. 
A comparison of this curve to the data has x 2 '/dof = 0.47 with 5 degrees of freedom. 



Compared to y 2 /dof = 0.68 with 8 degrees of freedom for QED alone, we conclude 
that this peak is not statistically significant. 

We extract 90% confidence level upper limits on (2 J + l)r ee BR 77 , for the nor- 
malization of 0.955 and particular choices of M and T, by numerical integration of 
the likelihood function.^ These limits are shown in Figure 2. It is apparent that the 
limits are insensitive to the total width T when it is comparable to or smaller than 
the center-of-mass energy spread a E . Furthermore, the limits shift only slightly if the 
normalization is varied from its maximum likelihood value. 

Upper limits on (2J+ l)r ee BR 77 at the 90% confidence level range from 0.5 keV 
to 8 keV, for an X mass between 57 and 59.6 GeV and a total width T below 100 MeV. 
Our limits can be compared directly with the expectation from the L3 observation 
of e + e~ — > £ + £~77 under the assumption that their clustered events arose from the 
production of a new state X. For example, if we assume that X is a scalar that 
couples only to photon and lepton pairs with a total width of T = T 77 + 3r ee , then 
an upper limit on (2 J + l)r ee BR 77 can be converted into an exclusion region in the 
plane of T 77 vs. r ee . Our 90% C.L. upper limit is shown in Figure 3 for a scalar X 
of mass 59 GeV. 

Our limits may be incorporated with 90% C.L. limits of other experiments such as 
the L3 measurement' 11 of e + e~ — > £ + £~ , yy and the OPAL measurement' 161 of e + e~ — > 
777. For the former case, we calculate' 141 the cross section using the automatic 

[15] 

amplitude generator GRACE including all possible diagrams. The limits from the 
three events with M 77 clustered near 59 GeV observed by the L3 collaboration are 

f Since the data are consistent with pure QED, the confidence interval includes the point 
(2J+ l)r co BR 77 = 0. 
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shown as the dashed curve in Figure 3, while the limits imposed by the absence of 
a invariant mass near 59 GeV are given by the dash-dotted curve. The limits 
set by the non-observation of an excess of events with M 77 = 59 GeV in the OPAL 
777 measurement are shown as the dotted curve of Figure 3. The limits of our 
measurement are thereby restricted to satisfying the condition that T 77 < 10.6 MeV 
and r ee < 120 MeV. To summarize, the four inclusion regions overlap in the lower 
left of Figure 3. 

For a (pseudo-) scalar that couples directly to pairs of fermions as well as photons, 
the exact forms of the decay widths for X — > e + e _ and X — > 77 are 

Q 2 

r ee = (6) 

8tt y ' 

and 

a 2 M 3 

F77 = 64^3 if (7) 
respectively, where g ee is the coupling constant of X to electrons, a is the fine structure 
constant, and Fx is a mass parameter analogous to f n in tt° decay. We use a /{it Fx) 
as the coupling constant of X to 77. Assuming that BR 77 ~ 1, we can convert the 
90% C.L. upper limit on r ee BR 77 from our likelihood analysis into an upper limit on 
g ee of 1.78 x 10~ 3 . We also place a 90% C.L. lower limit on F x of 730 MeV that comes 
from the OPAL upper limit on T 77 of 10.6 MeV. 1161 We examine the effect of these 
limits upon the eej vertex correction value, (g — 2) e , which will have an additional 
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contribution from this model of the form 

m e ag ee A 

Ax ^Fx-J^ ln M (8) 
where A is the ultraviolet cutoff (fixed at 1 TeV). Our limits on Fx and g ee im- 
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ply an upper bound of Ax < 2.0 x 10 , which is within the maximum allowed 
contribution [18] to (g - 2) e of 2.7 x lO" 10 . 

6. Conclusion 

We have searched for resonant production of a new state X in e + e~ collisions that 
couples to photon pairs at center-of-mass energies between 57.4 and 59.5 GeV. We 
find that the observed data are consistent with the QED prediction for e + e~ — > 77 
(x 2 /dof = 0.68 for 8 degrees of freedom). Furthermore, we have extracted 90% C.L. 
upper limits on (2 J + l)r ee BR 77 of 0.5 — 8.0keV for the process e + e~ — > X — > 77, 
should a state X exist with a mass between 57 and 59.6 GeV. 
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TABLE CAPTIONS 



Table 1. Results of the energy scan, y/s is the center-of-mass energy, a E is the 
width of the center-of-mass energy distribution, L is the integrated luminosity from 
the ESC, Nqed is the number of events expected from QED alone, iV 77 is the number 
of e + e~ — > 77 seen, and Ng B /L is the number of barrel Bhabha events seen per pb _1 
of luminosity, scaled by s. 



(GeV) 


a E (GeV) 


L (pb- 1 ) 


Nqed 


iV 77 


N^IN QED 


N BB /L 


57.374 


0.095 


2.05 


72.9 ± 0.4 


69 


0.95 ± 0.11 


265 ±11 


57.772 


0.096 


19.55 


680.5 ±3.4 


636 


0.935 ±0.037 


276.0 ±3.8 


57.972 


0.096 


1.39 


48.4 ±0.2 


49 


1.01 ±0.15 


272 ± 14 


58.220 


0.098 


1.33 


45.8 ±0.2 


50 


1.09 ±0.15 


259 ± 14 


58.470 


0.098 


1.33 


45.6 ±0.2 


54 


1.18 ±0.16 


308 ± 15 


58.718 


0.099 


1.73 


58.7 ±0.3 


61 


1.04 ±0.13 


278 ± 13 


58.968 


0.100 


1.36 


45.9 ±0.2 


45 


0.98 ±0.15 


284 ± 15 


59.216 


0.100 


1.21 


40.5 ±0.2 


45 


1.11 ±0.17 


294 ± 16 


59.466 


0.102 


0.97 


32.3 ±0.2 


26 


0.81 ±0.16 


236 ± 16 



FIGURE CAPTIONS 



Figure 1. Ratio of observed 77 events to the QED prediction. The solid curve is the 
QED prediction, with x 2 /dof = 0.68 for 8 degrees of freedom and a normalization of 
0.961. The dashed curve is the maximum likelihood prediction, with x 2 /dof = 0.47 
for 5 degrees of freedom and a normalization of 0.955. 

Figure 2. 90% C.L. upper limits on (2J+ l)r ee BR 77 as a function of mass M of the 
state X, for widths of 1 MeV (solid curve), 100 MeV (dashed curve), 1 GeV (dotted 
curve) and 3 GeV (dash-dotted curve), with a normalization of 0.955. The region 
above a given curve is excluded. 

Figure 3. Exclusion regions (90% C.L.) on T 77 vs. T ee for a scalar X with M ~ 
59 GeV, assuming T = T 77 + 3r ee - The region above and to the right of a given 
curve is excluded. The solid curve is from this experiment. The other curves are from 
the observation of up to three £+£"77 events with M 77 ~ 59 GeV (dashed) and the 
absence of £+£"77 events with Mi+£- ~ 59 GeV (dash-dotted) (Ref. 1), and from the 
lack of a signal at M 77 ~ 59 GeV in 777 events (dotted) (Ref. 16). 



